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Abstract 

Thermalization in an expanding parton plasma is studied within the frame- 
work of Boltzmann equation in the absence of any mean fields. In particular, 
we study the time-dependence of the relaxation time to the lowest order in 
finite temperature QCD and how such time-dependence affects the thermaliza- 
tion of an expanding parton plasma. Because of Debye screening and Landau 
damping at finite temperature, the relaxation time (or transport rates) is free 
of infrared divergencies in both longitudinal and transverse interactions. The 
resultant relaxation time decreases with time in an expanding plasma like l/r^, 
with P < 1. We prove in this case that thermal equilibrium will eventually be 
established given a long life-time of the system. However, a fixed momentum 
cut-off in the calculation of the relaxation time gives rise to a much stronger 
time dependence which will slow down thermal equilibrium. It is also demon- 
strated that the "memory effect" of the initial condition affects the approach 
to thermal equilibrium and the final entropy production. 
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1 Introduction 



Perturbative-QCD-based models developed in the last few years predict that nucleus- 
nucleus collisions at future collider energies are dominated by hard or semihard pro- 
cesses in, H, H, H . These processes happen during the very early stage of the collisions 
and they produce a rather large number of semihard partons which essentially form 
a hot and undersaturated parton gas @, ^. However, this parton gas is initially far 
away from thermal and chemical equilibrium [^]. Secondary parton scatterings in 
the gas may eventually lead to local thermal and chemical equilibrium if the parton 
interactions are sufficiently strong. 

Transport calculations based on a semiclassical parton cascade model indicate 
that thermal equilibrium could be established within a rather short time of about 
1 fm/c. However, the complexity of the Monte Carlo simulations makes it difficult 
to obtain a lucid understanding of the dependence of the thermalization time on the 
many parameters employed in the model. One such parameter is the cut-off of mo- 
mentum transfer in binary parton scatterings. The cut-off was first introduced to 
regularize the infrared divergency of the cross section between two massless partons 
in high-energy pp and pp collisions The value of the momentum cut-off is deter- 
mined phenomenologically to reproduce the measured total cross sections of pp and 
pp collisions. However, this cut-off is not necessary anymore in a high-temperature 
quark-gluon plasma, since the Debye screening and Landau damping provide natural 
regularizations of the infrared divergency. Since transport times depend sensitively 
on the screening masses which in turn depend on the temperature, the introduc- 
tion of an artificial cut-off could give rise to a completely different behavior of the 
thermalization time and consequently the approach to thermal equilibrium. 

The approach to thermal equilibrium in relativistic heavy ion collisions is dictated 
by the competition between expansion and parton interactions . If the expansion is 
much rapid than the typical collision time among partons, e.g., shortly after partons 
are initially produced, the expansion is closer to free-streaming than hydrodynamic 
expansion. Only at times in the order of the collision time may the parton gas 
reach local thermal equilibrium and expand hydrodynamically. Furthermore, the 
time dependence of the collision time (or the relaxation time) will determine whether 
the system can eventually reach local thermal equilibrium because of the competition 
between expansion and parton interactions. If the collision time increases rapidly with 
time, the parton system may never thermalize, leading only to a free- streaming limit. 
The collision time, therefore, is a very important quantity which in turn depends 
sensitively on the infrared behavior of parton interactions. 

In QCD, parton scattering cross sections exhibit a quadratic infrared singularity 
due to the exchange of a massless gluon. The infrared behavior can be improved by in- 
cluding corrections from hard thermal loops to the gluon propagators. Resummation 
of these thermal loops gives rise to an effective gluon propagator which screens long 



range interactions (Debye screening). Braaten and Pisarski have developed this 
resummation technique systematically and used it to calculate the damping rate of a 
soft gluon {p ~ gT) which is gauge invariant and complete to the leading order in the 
QCD coupling constant g [|lOl- For a fast particle (p^T), the exchanged gluons probe 



the static limit of the magnetic interactions which by the transversality condition are 
not screened. One thus has to introduce a nonperturbative magnetic screening mass 
to regularize the logarithmic infrared singularity in the static limit |Tl|, |12|, |l3l. As a 



result, the damping rate for an energetic particle to the leading order in g is 

T r^T[as\n{l/as) + 0{as)] , (1) 

where as = g'^/4:7i. However, as we will argue, damping rates do not determine how 
fast a system approach local thermal equilibrium. What really determine the thermal- 
ization processes are the transport rates which are free of the logarithmic divergency 
after the resummation of thermal loops [|1^, |T^. This is because thermalization is 



achieved to the leading order mainly through momentum changes in elastic scatter- 
ings. Thus, the effective cross section should be weighted by the momentum transfer 
and the dynamic screening due to the Landau damping of the gluons is sufficient to 
regularize the logarithmic singularity in the transverse interactions |T^, ^7\, [1^ . The 



resultant transport times for a system near thermal equilibrium behave like 

— ^TalHl/as) (2) 

to the leading order in as. 

For a system near local thermal equilibrium, the time dependence of the transport 
times is through the temperature according to Eq. (^. This dependence is in general 
slower than 1/r and thus can lead to local thermal equilibrium according to our 
earlier argument based on the relaxation time approximation [Q. However, if one 
introduces an artificial cut-off for the momentum transfers of elastic parton scatterings 
as in the numerical simulation of a classical parton cascade @], the time dependence 
will be much stronger. Consequently, as we will demonstrate in this paper, the 
system will approach local thermal equilibrium much slowly. We will also demonstrate 
that inclusion of the screening effects is the key to a slower time dependence of the 
relaxation time, therefore a faster approach to thermal equilibrium. 

This paper is organized as follows. In Section II, we first re-examine the Boltzmann 
equation and the evaluation of the damping rate and the relaxation time to the lowest 
order, including only 2^2 processes. We will also discuss the time dependence of 
the relaxation time in different scenarios. In Section HI, we will solve the Boltzmann 
equation in the relaxation time approximation and demonstrate how time dependence 
of the relaxation time will affect the approach to thermal equilibrium. We also show 
how initial conditions of a system affect the thermalization processes and the final 
total entropy production (or "memory effect") in Section IV. Finally in Section V 



we give a summary and an outlook, especially of the numerical simulations of parton 
thermalization, taking into account of the Debye screening and Landau damping 
effects without double counting. 



2 Time dependence of thermalization time 

In a system with two- components as, e.g., the quark and gluon plasma, the one 
that interacts the strongest will thermalize faster than the other. Subsequently there 
will be momentum and energy transfer between them. Since numerical simulations 
indicate that the initially produced partons are mostly gluon, we consider here a 
gluon gas only for simplicity. The Debye screening of color fields in the presence 
of semihard gluons |I9| will also allow us to neglect the effect of mean fields. We 
furthermore assume that the spatial variation of the system is small on the scale of 
a collision length so that we can approximate the evolution of the system by the 



Boltzmann equation |^ 



Vl-dfM = -Z/2/c^MMpA34[/]^|Mi2^34|'(27r)V(Pi + P2-P3-P4),(3) 
Pl234[/] = /l/2(l±/3)(l±/4)-/3/4(l±/l)(l±/2), (4) 

where Pi = {\pi\,Pi) are the four-momenta of massless partons and dpi = d'^p/{27T)^. 
To keep the formula general, ± are used for bosons (gluons) and fermions (quarks 
and anti-quarks), respectively. The statistical factor z/2 is 2{N^ — 1) = 16 for gluons 
and 12Nf for Nf flavors of quarks and anti-quarks with Nc = 3 colors. The squared 
matrix element, |Mi2^34p = 1/^12^341"^ / {I6E1E2E3E4) is summed over final states 
and averaged over initial states. For gluon-gluon scatterings, 

|-Mi2^34r = C^..4,^(3-f -5-^), (5) 

where Cgg = N^/{N^ — 1) = 9/8 is the color factor of gluon-gluon scatterings, s, 
t and u are the Mandelstam variables. There is clearly a quadratic singularity for 
small energy uj and momentum q transfers because of the long range interactions 
mediated by the massless gauge bosons. Because the final state has two identical 
particles, su/t"^ should contribute equally as st/u^ in Eq. (|^). We thus can approxi- 
mate |A1i2^34p ~ 'iCggQ^s'^ /t^ for small-angle gluon scatterings. Since the collisional 
integral is dominated by contributions from near the singularity, we can assume a 
small angle scattering approximation, i.e., uj,q ^ Ei,E2. Then energy-momentum 
conservation leads to 

P3 = Pi + q, P4 = P2 - q , 

E3 = El -\- LU, i?4 = E2 — UJ , 

LU ^ Vi-q^V2-q. (6) 



The integration over ^3 and ^4 can be rewritten as 



(27r)^y' dpsdp^6\Pi + P2-P3-P4) = J d?q J'^duj6{uj -Vrci)6{uj -v^-q) . (7) 

In a medium, one can use a resummation technique to include an infinite number 
of loop corrections to the gluon exchange. Using Dyson's equation, this amounts to 
an effective gluon propagator. One can use this effective propagator to obtain the 
effective matrix element squared for forward gluon scatterings (see Appendix A), 



\M, 



'1 — X ) COS( 



where cos0 = (vi x q)-(vi x q) and x = oj/q. The scaled self-energies in the long 
wavelength limit are given by | pT[ | 

1 + X' 
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(9) 
(10) 



where = g'^{Nc + A^//2)T^/3 is the Debye screening mass in thermal QCD. The 
imaginary parts provide Landau damping to parton interactions in a thermal medium. 
We see that the longitudinal interactions are screened by thermal interactions. How- 
ever, the transverse interactions still have a logarithmic singularity in the static limit. 
This singularity can only be regularized by introducing a nonperturbative magnetic 
screening mass in the calculation of the damping rate of a fast parton. 
We can use the definition of parton interaction rates. 



rioss(Pi) 

rgain(Pl) 



(27r)5 

J^2 

(27r)5 



J d^p2d'qduf2{l ± /3)(1 ± f4)^\Mu^u\'S{u - vi-q)5(^ - V2-®1) 
J d'p2d\duhUil±f2)^\Mu^u\^6iuj-Y^-q)6iuj-V2-q} , (12) 



and rewrite the Boltzmann equation as 
vi-dfi = 



where 



= -/irioss(j3l) + (1 ± /l)rgain(j3l) 

= -r(pi)(/i-/), 

T{p) = Tiossip) =F rgain(p) 



(13) 



(14) 



is usually referred to as the damping rate of a particle (or a quasiparticle), and / is 
defined as 

rgain(p) 



f{p) 



rioss(p) =F rgain(p) ' 



(15) 



For a system in local thermal equilibrium, one can relate the damping rate to the 
imaginary part of the gluon self-energy p2| , 



ImWJp) = -2{p-u)Tip) , (16) 



where the factor 2 comes from our definition of the interaction rates among identical 
particles. One can also show that, if f{p) takes the local equilibrium form f^^ip) = 
{exp{p-u/T) =F 1)"S 



rgain (p) 

using the energy and momentum conservation and the identity l±f'^'^{p) = /^^(p) exp(p 
u/T). Therefore, by definition, f{j>) becomes f^'^ij))- Thus, the global equilibrium 
distribution f^\p) is a solution to the Boltzmann equation if the flow velocity u is 
independent of space and time. 

We can complete the angular integrations in Eqs. (|ll]) and ([T^) . Making approx- 
imations /(ps) ~ /(Pi) and fipi) ~ f{P2), we can also complete the integration over 
P2 by cutting off the integration over q at gmax ~ 3T ^7^. We then obtain the gluon 
damping rate, 

including only gluon-gluon scatterings. 

The contribution from longitudinal interactions is finite and proportional to g^T 
due to the Debye screening. However, Debye screening is absent in the transverse 
interactions in the static limit. There is a logarithmic divergency even if Landau 
damping is taken into account. One solution to this problem is to add a nonpertur- 
bative magnetic mass rrimag ~ g^T to the transverse self-energy TCTix). In this case, 
the dominant contribution of the transverse interactions comes from m^i^g^q^qD 
and thus is proportional to g'^T\n{qD/m^s.g) which is independent of gmax ^ Qd- 
With TiL^x) and nT^x) given by Eqs. (|) and (p!OD, we can complete the numerical 
integration. A fit to the numerical result gives us 



■-gg ^ ^gg^g 

47r 6 qj^ 



(In ^^-1.0 + 2.0^^ -0.32^) + 1.1 



^mag to ?max Q'max ~^ 1d 



(19) 



where the first term comes from the transverse interactions while the second from the 



longitudinal ones. Using the estimate of mmag ~ 0.255'^ Nc/2g^T from Ref. |2^ and 
neglecting the quark contribution to the Debye screening mass, we have 

^ N,asT [ln(l/«,) - 0.1 + 0(a,)] . (20) 

Note that contributions to the order in Eq. (|l^) have been neglected, since they 
are not complete in our calculation. In order to have a complete calculation of such 



higher order corrections, one has to include thermal vertex and vacuum corrections 
which should depend on the renormalization scale. The final result to this order 
should be invariant under the renormalization group. This result agrees with previous 



calculations |Tl|, |12[ ^ to the leading order of which depends only on the imaginary 
part of the transverse self-energy, nxlx) ~ — ■i(7r/4)g|)X at small x. Inclusion of the 
full expression of the self-energy only contribute to the next order corrections. The 
increase in scattering by including quarks is exactly compensated by the increase in 
Debye screening due to quarks |]12| . 

For a soft gluon {p ~ gT), one can not neglect its thermal mass anymore. The 
damping rate for a gluon at rest will not have the logarithmic divergency in the 
transverse interaction, since the exchanged gluon must carry nonzero momentum and 
energy at least of order of gT and thus never approach to the static limit. In addition 
the transverse (magnetic) interactions are reduced by velocity factors which for the 
massive partons are smaller than the speed of light. The damping rate in this case 
was found by Braaten and Pisarski [|10| to the leading order as 

rf^'(O) ^ l.lNrasT . (21) 



Apparently, the damping rate has a nontrivial momentum dependence 

As we have mentioned, f^^{p-u/T) is a solution to the Boltzmann equation as far 
as the fiow velocity is uniform in space and time. We should emphasize here that the 
damping rate does not determine how rapidly a system near equilibrium approaches 
it as one would naively think. The thermalization time is actually related to the 
transport rates [|14|, |l8l. The easiest way to prove this is to check that 



the logarithmic divergency that has plagued the calculation of the damping rate of a 
fast gluon does not appear in the Boltzmann equation. To check this, we make the 
following expansion: 



f{P3A) ^ /(Pl,2) ± CU/'(P1,2) + y/"(Pi,2) , (22) 

for small angle scatterings. The function -Fi234[/] in Boltzmann equation becomes 

i^l234[/] = [/l(l ± fl)f^ + /2(1 ± /2)/r - 2(1 ± /i ± , (23) 

which is proportional to q^. Here we have dropped terms linear in x since they vanish 
after integration over x. One can verify that this function after expansion still vanishes 
for the equilibrium distribution, -Fi234[/'^^] = 0. For a system away from equilibrium, 
the coUisional integral in the Boltzmann equation is nonzero but finite despite the 
logarithmic singularity in the transverse part of the matrix element squared, |Mi2^34p 
when q,uj ^ 0. Because of the factor in -Fi234[/], Landau damping in the self-energy 
of the exchanged gluon is sufficient to give a finite value of the collisional integral. In 
other words, thermalization not only depends on the parton interaction rates but also 



on the efficiency of transferring momentum in each interaction. Those interactions 
with zero energy and momentum transfers do not contribute to the thermahzation 
process, though their cross sections are infinitely large. This is why the thermahzation 
time and other transport coefficients do not suffer from the infrared divergency as 



pointed out in a number of papers |12, 16, 17, 18 



For a system near equilibrium, one can characterize the deviation from equilibrium 
by 5/ = / — /'^^ = —6{p)p-df/{p-u) in a relaxation time approximation. Equivalently, 
one has 

—P-9f = --^- (24) 

In general the relaxation time 9{p) depends on momentum p and in principle can 
be obtained by solving the linearized Boltzmann equation. In this paper, we neglect 
the momentum dependence of the relaxation time. For a pure gluonic gas near local 
thermal equilibrium where thermahzation is achieved through viscous relaxation, the 
relaxation time is (Appendix B) 

1.0.92iVMl„(-H.). (25) 

Again, because of the extra factor appearing in the transport rate, the Debye 
screening and Landau damping are sufficient to regularize the effective transport 
cross section. The dominant contribution comes from interactions with q'D~Q'~Q'max, 
leading to a logarithmic factor ln(gmax/5'z)) as compared to ln(g£)/mmag) in the gluon 
damping rate. Therefore, the dependence of the relaxation time on the (weak) cou- 
pling constant and the color dimension is quite different from the gluon damping 
rate. 

If the system is close to thermal equilibrium, the hydrodynamic equations from 
energy-momentum conservation to the zeroth order of 5f can give us the time evolu- 
tion of the temperature T. For an ideal gluon gas with one- dimensional expansion, T 
decreases like T/Tq = {tq/tY/^. Therefore, the relaxation time 9 increases with time 
with a power of 1/3, 

0.92N^Toa^,ln{1.6/N,as) ' ^ ' 

A more general time dependence of the relaxation time can have a power-law form, 

= e^{T/Tof , (27) 

which also covers both the constant {(3 = 0) and the linear {f3 = 1) CcLSGS cLS have been 



studied by several authors ^ |2^. The latter case arises when a constant 

scattering cross section, a, is assumed for the relaxation time, 6 ~ {an)^^, and with 
a density decreasing as n ~ r^^ due to one- dimensional expansion. For a system 
far away from thermal equilibrium, the time dependence may differ from Eq. (pG]). 
Our earlier calculations pi show that an initially free- streaming system has only 



a logarithmic time dependence. This is because the phase space for small angle 
scatterings (g ~ qo) opens up quadratically with time in the free-streaming case and 
it balances the decrease in parton density. 

We would like to emphasize that the weak time dependence of the relaxation time 
in Eq. (|26D depends very sensitively on the Debye screening of the small-angle parton 
scatterings which restricts the momentum transfer to q^qo = qT. Smaller Debye 
screening mass due to the decrease of the temperature, gives a larger interaction rate 
which then compensates the decrease of the parton density and thus gives the weak 
time dependence of the relaxation time. If we use a fixed momentum cut-off gcut, as 
in most of the numerical simulations of parton production and cascade ^ , instead 
of a time dependent Debye screening mass, the effective transport cross section will 
remain constant, proportional to A^^a^/T^ ln(T^/g^^J. The resultant relaxation time 

-^^NlalT\n{T/q,^,) (28) 

will increase more rapidly with time. 

If we have to include transverse expansion later in the evolution of a system, 
then the temperature will decrease faster, like tq/t, than in the one- dimensional case 
assuming hydrodynamic expansion. The relaxation time even with the inclusion of 
Debye screening will increase linearly with time. The relaxation time with a momen- 
tum cut-off applied to parton interactions will increase faster than linear with time, 
which will only lead the system into free-streaming. 

At this point we should emphasize that we have only considered the lowest or- 
der contribution from 2^2 processes in our calculation of the relaxation time. In 
principle, higher order processes, like 2 ^ 2 + n, should also contribute to the ther- 
malization. Such processes can be included by considering high order thermal vertex 
corrections. For a complete calculation, one should also include vacuum corrections 
and the result should depend on the renormalization scale and obey the renormaliza- 
tion group equation. In general, contributions from 2 ^ 2 + n processes should have 



a form 29 



^„,~^o[a,ln(gVg2)]^ (29) 

where q^ is the momentum scale of these processes. At zero temperature, go is some 
confinement scale below which perturbative QCD is no longer applicable. At finite 
temperature, go is very likely to be replaced by screening masses. Since the largest 
momentum scale in a system at finite temperature is g^ ~ T^, the leading correction 
from 2 ^ 2 + n processes must be, 

r„~ro[a.ln(l/a,)r. (30) 



Such corrections therefore are high orders in ln(l/Q;s) and are negligible in the week 
couple limit. For temperatures not far above the QCD phase transition temperature 



Tc ~ 200 MeV, the strong coupling constant is not very small. The above contribu- 
tions might not be negligible. However, for an order-of-magnitude estimate, we can 
neglect these higher order contributions. If one considers the chemical equilibration 
of a kinetically thermalized system as in Ref. n ^ m multiplication processes 
become very important. The leading contribution to the chemical equilibration, in 
this case, comes from 2^3 processes. 



3 Approach to thermal equihbrium 

Let us consider the early stage of a very heavy ion collision where transverse expansion 
is not important yet. We then can treat the system as a one-dimensional system. We 
assume along with Bjorken a scaling flow velocity 



X 



Un. = — = (coshr], sinh?7, 0^) (31) 
r 

in the longitudinal direction, where 

r = V^^, r^=^\n(i±^^ (32) 

are the proper time and spatial rapidity, respectively. In terms of these new vari- 
ables, the Boltzmann equation in the relaxation time approximation for a system 
near thermal equilibrium becomes 

df tanhe^/ f-r 



where = 1] — y and y is the rapidity of a particle. 



(33) 



^-^InlM. (34) 




Since p-u = prcosh^, we can see that the solution to the above Boltzmann equation 
is a function of ^, r and px, and is Lorentz invariant under longitudinal boost. One 
may then, as done in |^, solve the Boltzmann equation in the central slice only, 
i.e., rj = 0. To simplify the problem, we assume the system is already in chemical 
equilibrium so that the gluon chemical potential vanishes. This is not always true in a 
realistic situation as shown by numerical simulations of initial parton production 
at around RHIC energies. However, at LHC energies, the small-x behavior of the 
parton distributions as measured by recent HERA experiments gives much higher 



densities of initially produced partons very close to chemical equilibrium For 
a given momentum and time dependence of the relaxation time, one can find the 
solution to Eq. (^3]) in an integral form. 



f{PT,^,r) = e>^»->^/o(pT,sinh-i(^^)) 

+ rdx'e^-^'r{pT,smh-\^^^^),T')- (35) 



Here /o(pt,0 is the initial distribution at time tq, and the time dependence of the 
temperature in f^^pr-, ^, T) is determined by requiring the energy density for / and 
f^"^ be equal at any time, 

e(T) ^ J dpEfip) = J dpEr{p,T) . (36) 

The variable x is defined by 

x(PT,e,r) = Xo+ rrfr'ri(pT,sinh-i(^^),r'). (37) 

J TO T 



If we take the form in Eq. (|27D for a momentum-independent relaxation time, then 

X = Xo (—) , Xo = , (38 

Vro/ 1-/3 



for /3 7^ 1. The case for P = 1 was studied by Gavin p8|. For a more general 
discussion, let us also assume a initial parton distribution 

= (39) 

with a simple plateau distribution in ^. One can also consider a gaussian form of 



distribution |]3^, but the above form is simpler. The width F is a measure of the 
initial correlation between space and momentum, or rj and y. This width which 
normally depends on the transverse momentum px 0] can be estimated by the 
uncertainty principle. We will assume it a constant for simplicity and will study 
how the thermalization process depends on the initial condition by varying Y. Here, 
w{pt) is the pt distribution which usually has a power-law or exponential form. The 
corresponding phase space density within a volume V is 

MpTX) = '-^^%m\-Y), (40) 

Ug V Pt cosh 4 

with an initial energy density 

eo = ^ ^^"^ ^ J dplpTwipr) ■ (41) 

Following Baym we study the solution to the Boltzmann equation by taking 
the first moment (energy density) of Eq. ( P^ ) with respect to the single parton energy. 
Defining, 



re r 



G = ^ , (42) 
Toeo 



Gix) = e^'^'^Hyiixo/x) — ) + / dx'e^ "^G(x')/^((x7x)~) , (43) 



we have from Eq. (^ 



where x'/x — i'^' I'^Y ^ from Eq. (p^), and 

Hy{a) = — ^ di^l + aHin\i^i, (44) 
, / \ 1 / arcsin Vl - 

If the initial distribution /o is an equihbrium one, one should replace -ffy(a) by h{a) 
in Eq. (^). However, Hyia) [with -ffy(l) = 1 and -ffy(O) = y/sinhy] covers a 
wide selection of initial distributions with different values of Y. For Y = 1.23 (i.e., 
Y/ sinhy = h{0) = vr/4), ify(a) is very similar to h{a) for an isotropic initial dis- 
tribution. When Y = 0, the initial distribution corresponds to the Bjorken scaling 
ansatz ij] = y). One then has Hq^o) = 1, which is similar to the case discussed by 
Baym |^. After performing a partial integration in Eq. (^), we have finally 

rdx'e^'-f^[Gix')hi{x/x)^)] = e^^{Hy{T,/r) - h{r,/r)] . (46) 



The discussion of the condition for equilibrium at large times, r, is now similar to our 
earlier analysis 0] even with the general initial condition assumed. However, as we 
will show later, the initial condition will have important influences on the approach 
to equilibrium and the final entropy production. 

For /? < 1, r ^ oo corresponds to x oo. In this limit, the r.h.s. of Eq. ( ^Bj ) is 
a finite number, e^°{Y/ sinhy — vr/4). Since the integrand on the l.h.s. of Eq. (^6|) 
has an exponential factor, a finite integral must require 

^[G'(x')M(x7x)^)]x'=x = 0, X-oo. (47) 

Using h{l) = 1 and h'{l) = 1/3, one can find the solution to the above equation, 

G (X oc (ro/r)^/3 , t/tq ^ oo , (48) 

which corresponds to the one- dimensional hydrodynamic limit, e oc (tq/t)^/'^. There- 
fore, thermal equilibrium will eventually be established for j3 < 1. 

For /3 > 1, r ^ oo limit corresponds to x —0. Using h{0) = 7r/4, one has from 
Eq. (ID 

Gix) = e^" - 7 r dx'Gix') = const. , x - -0 • (49) 
4 JO 

This corresponds to the free-streaming limit, e oc tq/t. Therefore, thermal equilib- 
rium will never be achieved if /5 > 1. 



The special case /? = 1 was studied by Gavin |^ who found that the system will 
also reach an asymptotic state lying in between free-streaming and hydrodynamic 
limit depending on the value of the prefactor 6q in Eq. (|27|). In the asymptotic state. 



G{t/tq) oc (tq/t)^ with < 7 < 1/3. Only for a very large 9o/tq ^ 1, does the 




Fig. 1 

Figure 1: Time evolution of G{t /tq){t /tq)^/^ = (e/eo)(r/ro)'^/'^ according to the so- 
lution to Boltzmann equation, for different time dependence of the relaxation time 
e = Oq^t/tqY, with eo/TQ = 1, /5 = 0.9 (dotted), 1/3 (solid), 1/4 (dot-dashed) and 
(dashed line). Yq = 1.2 is the width of the initial rapidity distribution. 

system approach to the hydrodynamic limit with 7 ^ 1/3 ~ 16ro/45^0) which is very 
close in form to the case of /3^1 as noted by Gavin in Ref. [p^ . 

Shown in Figs. |l| and |^ are the numerical solutions of Eq. (^3]) for different re- 
laxation times. For a clear presentation, we plot [r /tqY/'^G{t /tq) = (r/ro)'^/'^(e/eo) 
as a function of t/tq for different values of (3 and Oq/tq. In these two plots, we 
have chosen Y = 1.2 which corresponds to an isotropic initial distribution in mo- 
mentum space. If the system undergoes free streaming, i.e., G{r/To) = const., the 
plotted quantity should increase with t/tq with a power of 1/3. Thus we see from 
Figs. |l| and ^ that a system must undergo free streaming for a period of time before 
it approaches the equilibrium limit when parton interactions eventually balance the 
expansion. The duration of such a period, which we can define as the thermalization 
time Tth is determined by the relaxation time and its time dependence. For fixed (3, 
the thermalization time Xth is approximately proportional to ^^o- From Fig. ^ we can 
estimate that Tth — S^o for P = 1/3 if we consider that thermalization is reached 
when G{t/tq) is about 10% from its hydrodynamical limit. Assuming one scattering 
is sufficient to thermalize the system |^ is therefore a serious underestimate. When 



{3 is close to 1, the system approaches the hydrodynamical limit very slowly as seen 
in Fig. |l[ However, the hydrodynamical limit is still achieved as long as /5 < 1, unlike 
when (3 > 1. When the thermalization time is very long, one must also consider 
whether 3-dimensional expansion and/or hadronization occur earlier. 




Fig. 2 

Figure 2: Same as Fig. (0), expect for P = 1/3, Oq/tq = 2 (dashed), 1 (solid) and 0.5 
(dot-dashed hne). 

4 Memory effect 

Let us now discuss how the initial condition /o influences the thermalization process 
and the flnal entropy production. In Fig. ^ the time evolution of {t/toY^^G^t/tq) is 
shown for different values of Y which characterizes the initial distribution of partons 
in the phase space. For large values of Y, we notice that the system initially expands 
even faster than the ideal hydrodynamical case and then turns over, approaching the 
hydrodynamical limit. To understand this, let us take the flrst moment of Eq. (|33D 
with respect to the single parton energy. We then have, by energy and momentum 
conservation, 

^ + ^ = 0, (50) 

dr T 

where the energy density is deflned by Eq. (136|) and 

is the longitudinal pressure. The solution to the Boltzmann equation can be parametrized 
as PL(r) = 7(r)e(r) with < 7 < 1. We have then from Eqs. (|50|) and (||), 

G(r) = e J^o''^^ . (52) 



Thus, 7(r) characterizes the thermalization of the system (or equal partition in lon- 
gitudinal and transverse direction). In the hydrodynamical limit, 7 = 1/3, while 



1.6 




Fig. 3 

Figure 3: Same as Fig. ([l|), except for P = 1/3, Oq/tq = 1, Yq = 0.2 (dashed), 1.2 
(solid), 2.5 (dotted) and 5.0 (dot-dashed hne). 



free-streaming corresponds to 7 = 0. One can expect that the initial evolution of the 
system near t/tq ~ 1 should be determined by the value of 70 = 7(to), 

G(r) ^ (ro/r)^o , t/tq ~ 1 . (53) 

Using the initial distribution in Eq. (^0]), we have 



Piiro) arctan(sinhr) 
^° = ^ = '- sinhF ■ ^^'^ 

An isotropic situation corresponds to arctan(sinh Fq)/ sinhyg = 2/3 or Yq = 1.167 
which is very close to the value we obtained by requiring ify(O) = tt/A. For Y < Yq, 
7o < 1/3, the system starts its evolution more like free streaming as we have noticed 
in our numerical solutions. The extreme case is the Bjorken scaling ansatz, Y = 0, 
7o = 0, which corresponds exactly to free streaming. For Y > Yq, 70 > 1/3, the 
system will initially expand in the longitudinal direction even faster than a thermal 
expansion, as also demonstrated in our numerical solutions. Physically, this is caused 
by the higher longitudinal pressure built up by the large amount of partons which are 
distributed over a large range of rapidity in the local frame. One can check that the 
expansion in this early stage is still dominated by free-streaming. However, working 
against such high pressure costs energy thus leading to less entropy production as we 
now show in the following. 

Since the system will eventually approach the thermal equilibrium limit when 
/3 < 1, we can define a prefactor Aoo by 

G(r/ro)=^oo(ro/r)^/^ t/tq ^ 00 . (55) 
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Fig. 4 

Figure 4: The prefactor, = / / TqY^'^ {t~ / t~o oo), as a function of the Yq; the 
width of the initial rapidity distribution, for different relaxation times 9 = 9q{t/tqY 
(/? = 1/3), 6'o/to = 2 (dashed), 1 (solid), and 0.5 (dot-dashed line). 

This prefactor in general can only be calculated numerically and will depend on the 
relaxation time as well as the initial condition as we have seen in Fig. ^. Using the 
above expression, one can calculate the final entropy of the system per comoving 
volume P 

4ToeoTo 4 eo . , . 

TSf = -—^—G{t/to) = T0-—A00 ■ (56) 

6 iQ 1 6 iQ 

which is exactly proportional to • Here Tq is only a parameter in the time depen- 
dence of the final temperature, T = To(ro/r)^/^. Since the initial system is not in 
equilibrium especially for large values of Y, (4/3)eo/To should not be considered as 
the initial entropy density. Therefore, is not the absolute increase of the total 
entropy over its initial value, except when the initial distribution is an equilibrium 
one. Nevertheless, Aoo still carries a lot of information about the thermalization pro- 
cess of the system. We plot this prefactor in Fig. § as a function of Y for /3 = 1/3 
but for different values of 9o/tq. It is clear that there is less entropy production for 
larger values of Y, since the system has to work against increasingly high longitudi- 
nal pressure thus converting its kinetic energy to expansion energy. The system has 
the maximum entropy production when the initial distribution is that of the Bjorken 
scaling ansatz, Y = 0. In this case, one can calculate [§] that the entropy production 
increases with the relaxation time like oc {9o/ToY^'^^^~f^\ For very large values 
of Y, Aoo decreases shghtly with 9o/to as indicated in Fig. |[ since the system has 
to spend longer time work against the extraordinaryly high longitudinal pressure. It 
is the competition between long thermalization time (thus entropy production) and 



long duration of work against high pressure that reverses the 6q/tq dependence of 
entropy production ^oo- 



5 Conclusion and outlook 

In this paper, we investigated the thermahzation process in a one- dimensional expand- 
ing parton plasma within the framework of the Boltzmann equation. In particular, 
we have studied the time dependence of the relaxation time and its influence on the 
thermahzation. If the time dependence is weaker than a linear form, we find that the 
thermal equilibrium limit will eventually be reached. For a time dependence stronger 
than the linear one, the system will never thermalize, only leading to a free-stream 
limit. For an exact linear time dependence, the system will reach an asymptotic state 
between free-streaming and thermal equilibrium. We find that the thermahzation 
process also depends on the initial condition of the system. The deviation of the 
initial momentum distribution from an isotropic one in the longitudinal direction de- 
termines the initial approach to thermal equilibrium. This initial approach will then 
carry its inertia throughout the whole thermahzation process. This "memory effect" 
can be seen from the dependence of the final total entropy production on the initial 
momentum distributions. 

Using perturbative QCD at finite temperature in the transport theory, we have 
calculated the relaxation time as a result of parton scatterings. We have pointed out 
the important differences between parton damping rates and thermahzation times. 
To regularize the singular behavior of the parton scattering cross sections, we have 
used a full gluon propagator which includes the resummation of an infinite number of 
hot thermal loops. Because of the singular behavior of the parton scattering matrix 
elements, the resultant relaxation time depends sensitively on the Debye screening 
mass qo = gT. In an expanding parton gas, the temperature decreases with time 
and so does the Debye screening mass, thus leading to an increasing transport cross 
section. This then compensates the decrease of the parton density and gives us 
a relaxation time with a weak time dependence, 9 oc {t/toY^^. However, if one 
introduces a fixed momentum cut-off to parton scatterings as in most of numerical 
simulations, one will introduce an extra logarithmic time dependence which will slow 
down the thermahzation process. 

Although we have demonstrated the necessity of the inclusion of Debye screening 
and Landau damping in the study of the parton thermahzation in an expanding 
system, it is not clear to us how to incorporate them into numerical simulations such 
as parton cascade models. One can include the Debye screening semiclassically 



by brutal force. However, one immediate problem we have to solve is how to avoid 
double counting. 

Since we used the full propagator which includes many thermal loops, we have 



effectively included multiple particle scatterings. One can easily see this by expanding 
the full propagator in terms of the bare propagators at zero temperature and the self- 
energy from thermal loops. The contribution from the real part of the self-energy 
(mainly Debye screening) corresponds to multiple particle scatterings in the thermal 
bath. One can in principle include particle radiation and absorption by considering 
the thermal loop corrections to the full vertices. In doing so, one can automatically 
avoid both the infrared and coUinear divergencies which one normally regularizes 
by resorting to two additional cut-offs [p^] . 

As we have demonstrated in Appendix B, the imaginary part of the self-energy also 
contributes to the effective parton scattering in a parton gas. In fact, the imaginary 
part, which is responsible for Landau damping, is necessary to regularize the trans- 
verse interaction since there is no magnetic screening in QCD. The contribution from 
this imaginary part to an effective two-parton scattering corresponds to independent 
scatterings of the two partons off particles in the thermal bath. This can easily cause 
double counting in a parton cascade model because independent parton scatterings 
are also simulated over the volume of the system. One might be able to avoid this 
problem by introducing a length scale in the order of 1/gT, within which multiple 
particle scatterings are not allowed and only one effective two-parton scattering (with 
Debye screening and Landau damping) is permitted. 
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Appendices 

A Effective matrix elements of forward parton scat- 
terings 

In this appendix we separate the QCD interactions into longitudinal and transverse 
parts with inclusion of screening in these two parts and show the simplifications that 
appear in the limit of forward scatterings or small momentum transfer interactions. 
We adopt the same convention as in Ref. |^T| and split a vector into its components 



parallel and orthogonal to the flow velocity m^, {u^Uf^ = 1), such that 

uj = Q-u, Q^, = Q^,- u^,{Q-u) . (57) 



We can then denote a vector hj Q = [uj,q\, with Q'^ = u"^ — and = —q^. In 
the local frame, where the flow velocity is m = (1, 0), u and q become the time and 
spatial components of the vector. Similarly, one also defines a tensor orthogonal to 

u^,, 

g^iu = g^iu - u^,u^ . (58) 

The full gluon propagator with momentum Q is obtained from the vacuum polariza- 
tion by using Dyson's equation, 

where a is a gauge fixing parameter, the longitudinal V^" , and transverse tensor Vj^ 
are defined as 

= 7^w^-qv)m^-qv), (60) 

q^ 

which are orthogonal to and also to each other, i.e., 

g,.pr = Q,K" = K.v^T = 0. (62) 

In additon, one also has V^^Vpv = V^. The free gluon propagator at zero temperature 
is in this case 

D^'^{-,r + o9^)l,. (63) 
If we choose Feynman gauge (a = 1), the full propagator also satisfies Q^A'^'^ = 0. 



The transverse and longitudinal self-energies are |2l 

Ul{Q) = (1 - x^)7rL{x) , Ut{Q) = 7Tt{x) (64) 

with the scaled self-energies, 71l{x) and ttt{x), given in Eqs. (|^) and (p!OD, where 
X = oj/q. 

Using the full gluon propagator, A^*^, we can obtain the effective matrix element 
of quark scatterings qi{Pi) + qj{P2) ^ ^^(^s) + QjiPi) {i J), 

\Mgg\^ = C,ggH{4\P^■A■P,\'-2{P^■Ps){P2■A■A*■P,) 

-2{P2-P,){PvA-A*-P,) + {PrP,){P,-P^)\A\'} , (65) 

where Cqq = {N^ - l)/4:N^ = 2/9 is the color factor and 

A^/'A*'' = — \ — f66) 

|Ap = A^'^A* = + i . (67) 



For small angle scatterings, —Q'^/2 = Pi-P^ = P2-Pa ^ -Pi'-Pj = Ps-Pa, only the first 
term, corresponding to t-channel scattering, in Eq. (p5| ) is dominant. Furthermore, 
uj = xq ~ Vi ■ q ^ V2 ■ q from the energy and momentum conservation. In this 
approximation, one can verify that 



PvVl-P2 
Pi-Vt-P2 



Define cos i 



EiE2{l-x') . 
(vi X q)-(vi X q), we can express Vi-V2 as 



(68) 
(69) 



V1-V2 = (vi-q)(v2-q) 



cos I 



l-(vrq) 



lV2-qj 



X COS ( 



We have then 
\M 



C,,g^lQ{E^E2 



l-x^ 



X cos ( 



(f — Hi up- — (f- — IIt 



(70) 



(71) 



The matrix elements for gluon-quark and gluon-gluon scatterings are similar in the 
small angle approximation, except the color factors, Cgq = 1/2 and Cgg = N^/{N^ — 
1) = 9/8. For scatterings of identical particles, one should also multiply a factor of 2 
to take into account of the equal contributions of t and w-channel scatterings. Using 



Es ^ El, E4 ^ E2, and Eq. 



one arrives at Eq. 



B Calculation of the relaxation time 

For a system near thermal equilibrium, the energy-momentum tensor is given to the 
zeroth order of the deviation 6f by 

Tj^ = (e + P)n^«'^-(7^^P, (72) 

where e is the energy density and P the pressure. One can also split the derivative 
dfj, into components parallel and orthogonal to the flow velocity u'^, 

= u^D + 5^, (73) 

where D = u-d and d^ = d^ — u^D. The energy-momentum conservation dJI^^!^^ = 
can be rewritten as 

De+ {e + P)d-u = 0, 
(e + P)Du^ - d^'P = . (74) 

Conservation of entropy requires that the entropy density s = {e + P)/T fulfills 
d^{sUfji) = Ds + sd-u = 0. Thus we obtain 



where dP/de =1/3 for an ideal gluon gas. 

We assume the space-time variation of the distribution / is smalL To the leading 
order in a gradient expansion in the relaxation time approximation, we have 



5/ 



d{p) 



p-df 



eq 



9{p) 



r{i±r)p-d 



p-u 



(76) 



p-u p-u 

Using the hydrodynamical relations, Eqs. ([7^) and (|75D , one can rewrite the above as 



5/ 



(77) 



where U^y = {l/2){d^Uy + d^u^), = 9^P/(e + P) — d^T/T which vanishes for an 
ideal quark and gluon gas and S^i, is another tensor orthogonal to the flow velocity. 



To leading order in the gradient expansion, the Boltzmann equation becomes 

/r(i ± n^^-'ipi) = dp,dpsdp, [^(pi)sr + ^(P2)sr - e(P3)sr - ^(P4)sr] (79) 



/rvi'^(i ± /r)(i ± /r)-iMi2-.34p(27r)45^(Pi + P2-P3- p,) . 

In general, the relaxation time is momentum dependent and it, or the deviation from 
local equilibrium, must be determined by solving the Boltzman equation. As a first 
approximation, we assume 9 is independent of the momentum p (one can also assume 
9 (X p which is found to be a good approximation in calculations of viscosities [0). 
Multiplying both sides of the above equation by Y^^^{pi) and integrating over dpi, we 
have then 

!dpn{i±DT.\p) 



9 



0) 



where 

m 



(27r)8 



J d'pj^i ± /D / d'p^mi ± /D J d\ J du^ii:, + S2 - S3 - S4)^ 



[1 — X ) COS( 



6{uj — Vi-q)5(co' — V2-q) 



where again we have made the approximation of dominance of small angle scatter- 
ings in the local frame, in which u = (1, 0) and S^i, only has nonvanishing spatial 
components. 



= - [PiPj - -E Si 



and 



2 ^9 
3 



^2) 



Using Eq. (|^) we have 



(Si + S2 - S3 - S4)i, ^ iv2 - ViYq^ + {V2 - ViYq' - uj{vivi 



f if 



1^1; ) 



^3) 



for small q and uj. In terms of cos0 and x = uj/q, 

(Si + S2 - S3 - S4)^ = 2(1 - Vi-V2)[2g2 -2cjq-(vi + V2) +cj2(l + vi-va)] 

= 2g2(i-cos0)(l -x2)2(2-x2 + x2cos0) . (84) 

One can easily complete the angular integrations in £[S] similarly as in the calculation 
of the gluon damping rate, and has 

m = / dp^plfni ± /r) / dp,plfP{l ± /r)X(g„.ax, qo) , (85) 



(2vr) 

X(gmax,gD) = dq J ^dxq{l-x 



- 2 / 7 2/1 2\2 I 4**' 



|g2 + 7ri(a;)|2 



|g2(l — a;^) + 7rr(a;)p + 7r2(a;)] [g2(l — x"^) + 7r'r(x' 

After performing integrations like 



for a Bose- Einstein distribution f^^, where ({n) is the Riemann's (^-function, we obtain 
for the relaxation time, 

The integral X(gmax, (Jd) has no dimension and therefore should be only a function of 
gmax/^D- Numerical evaluation of the integral gives 

X(g„ax/gD) = 2.3 HqlJql) - 0.62 + 2.7^^ + 0(5^?^ (89) 



max 



for large values of gmax/ qo- Considering Debye screening from only gluon interactions 
we have 

1^ (^) + ^(«/) 



9 



Similarly as we argued after Eq. (pO|), contributions to the order must be neglected 
since they are not complete in our calculation. 
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